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Abstract 


In this paper, we study five luminous supernovae (LSNe) Ibe (SN 2009ca, ASASSN-15mj, SN 2019omd, 
SN 2002ued, and SN 2021bmf) whose peak absolute magnitudes Mpeak are ~ —19.5 to —21 mag by fitting their 
multi-band light curves (LCs) with different energy source models. We find that SN 2009ca might be powered by 
the °°Ni model since the required °°Ni mass (0.56 M.) is comparable to those of energetic SNe Ic, while the rest 
four SNe cannot be accounted for the *°Ni model since their derived ?9Ni masses are 21 Mo or the ratios of the 
S°Ni mass to the ejecta mass are larger than 0.2. This indicates that some LSNe might be powered by *°Ni decay, 
while most of them need additional energy sources. We then use the magnetar plus °°Ni model and the fallback 
plus °°Ni model to fit the LCs of the four LSNe that cannot be explained by the °°Ni model, finding that the two 
models can account for the four SNe, and the derived parameters are comparable to those of LSNe or 
superluminous SNe in the literature, if they were (mainly) powered by magnetars or fallback. We suggest that the 
magnetar plus °°Ni model is more reasonable than the fallback plus °°Ni model, since the validity of the fallback 
plus °°Ni model depends on the value of accretion efficiency (7) and favors a large 7 value, and the magnetar plus 
°Ni model yields smaller x? /dof values. It should be pointed out that, however, the fallback plus °Ni model is 
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still a promising model that can account for the four SNe in our sample as well as other LSNe. 
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1. Introduction 


Core-collapse supernovae (CCSNe) result from the explo- 
sions of massive stars with zero-age main-sequence mass 
(Mzams Z 8.0 Mo) (Woosley et al. 2002; Janka et al. 2007). It 
is believed that the progenitors losing most or all hydrogen 
envelopes can produce stripped-envelope SNe (SESNe) which 
can be divided into type IIb SNe, type Ib SNe, and type Ic SNe 
(see Filippenko 1997; Gal-Yam 2017 for reviews), for which 
the progenitors lost most hydrogen envelopes, all hydrogen 
envelopes, and most helium envelopes, respectively. 

Most types Ib and Ic (type Ibc) SNe are dimmer than type Ia 
SNe whose peak absolute magnitudes Mpeak are ~ —19.5 mag. 
Over the past two decades, however, a few hundred super- 
luminous type Ibc SNe (SLSNe Ibc) with Mpeak S —21 mag 
were confirmed and comprehensively studied (see Gal-Yam 
2012 and Gal-Yam 2019 for reviews). Additionally, a few 
dozen luminous type Ibc SNe (LSNe Ibc) whose Mpeak are 
between ~ —19.5 and ~ —21 mag are also been studied.? 


5 Richardson et al. (2002) show that Mpeg of SNe Ia are ~ —19.5 mag, 
while Gal-Yam (2012) suggests that the threshold of SLSNe is —21 mag. 
Therefore, Wang et al. (2015b) suggest that the SNe with Mpeak in —19.5 to 
—21 mag are LSNe. 


Unlike normal type Ibc which is believed to be (mainly) 
powered by the cascade decay of *°Ni synthesized by the 
explosions (Colgate & McKee 1969; Colgate et al. 1980; 
Arnett 1982), the energy sources of LSNe Ibc and SLSNe Ibc 
are still elusive. 

Previous studies for SLSNe Ibc indicate that almost all 
SLSNe cannot be explained by the °°Ni model (including the 
pair-instability SNe model which need a few M. of ?9Ni). To 
account for the light curves (LCs) of SLSNe which cannot be 
explained by the *°Ni model, three models which are the 
magnetar model (Maeda et al. 2007; Kasen & Bildsten 2010; 
Woosley 2010; Inserra et al. 2013; Wang et al. 2015a), the SN 
ejecta-circumstellar medium (CSM) interaction (CSI) model 
(Chevalier 1982; Chevalier & Fransson 1994; Chatzopoulos 
et al. 2012, 2013), and the fallback accretion model (hereafter 
the fallback model; Dexter & Kasen 2013; Moriya et al. 20182, 
2018b; Anderson et al. 2018; Li et al. 2020) have been 
proposed. 

The case for LSNe Ibc is more complicated. Wang et al. 
(015b) suggest that the magnetar plus °°Ni_ model is a 
promising model to account for LSNe Ic since the ?9Ni masses 
needed by the ?5Ni model are larger than the reasonable values, 
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Table 1 
The Information on Five Luminous Supernovae 
Name R.A. Decl. Type Redshift Meest.peak References? 
(72000) (72000) (mag) 
SN 2009ca 21°26™225 20 —40?2'48"6 Ic-BL 0.096 r= —19.77 1,2 
ASASSN-15mj 14^02"155 64 +33°2/40"29 Ib 0.034 V= —19.5 3 
SN 20190md 00^367575972 —03?2/39"17 Ib 0.082 r — —19.76 4 
SN 2021bmf 16^33"295 408 —06?2/49"58 Ic-BL 0.04 r = —20.4 4 
SN 2022ued 08"54™45 § 854 4-78?2/48"43 Ib 0.1087 r= —20.5 4 
Note. 


è References. (1) Taddia et al. (2018); (2) Stritzinger et al. (2018); (3) Shappee et al. (2014); (4) https: //sites.astro.caltech.edu/ztf/bts (Perley et al. 2020). 


while the contribution of a moderate amount of *°Ni cannot be 
neglected. Wang et al. (2015b) use the magnetar plus ?9Ni 
model to fit the bolometric LCs of three LSNe Ic-BL 
(SN 2010ay, SN 2006nx, and SN 14475), and find that the 
model can explain their LCs and the derived parameters for 
magnetar were Py ~ 7-15 ms and B ~ (4-7) x 10'* G. Gomez 
et al. (2021) investigated SN 2019stc and found that the first 
and second peaks in this SN could be powered by the magnetar 
plus ?9Ni decay and the ejecta-CSM interaction, respectively. 
Gomez et al. (2022) used the magnetar plus ?6Ni model to fit 
the multi-band LCs of 40 LSNe with r band Mpeak between 
—19mag and -—20 mag, finding that the parameters of 
magnetars were Po ~ 1-23 ms and B ~ (0.1-13) x 10!^ G. 

In this paper, we collect the multi-band data of five LSNe Ibc 
with Mpeak ~ —19.5 to —21mag, and constrain their physical 
properties by using the °°Ni model, the magnetar plus Ni 
model, as well as the fallback plus 56Ni model to fit their multi- 
band LCs. In Section 2, we model the multi-band LCs of the 
five LSNe using these three models. In Section 3, we discuss 
our results and draw some conclusions in Section 4. 
Throughout the paper, we assume €), = 0.315, Qa = 0.685, 
and Ho— 67.3 km s! Mpc | (Planck Collaboration et al. 
2014). The values of the Milky Way reddening (Eg. y) of all 
events are from Schlafly & Finkbeiner (2011). 


2. Modeling the Multi-band Light Curves of 
Five LSNe 


We collect about 50 LSNe Ibe with Mpeak ~ —19.5 to 
—21mag from the Open Supernova Catalog (Guillochon et al. 
2017)* and the public catalog of transients from the Zwicky 
Transient Facility (ZTF) Bright Transient Survey (BTS) (Perley 
et al. 2020)? We then exclude the LSNe which had been fitted 
by the literature (e.g., Wang et al. 2015b; Gomez et al. 2021, 
2022), and obtain the rest five LSNe, which are SN 2009ca, 
ASASSN-15mj, SN 20190md, SN 2021bmf, and SN 2022ued. 
Among the five LSNe, ASASSN-15mj, SN 20190md, and 


3 https: / /sne.space 


> https://sites.astro.caltech.edu/ztf/bts /bts.php 


SN 2022ued are SNe Ib, while SN 2009ca and SN 2021bmf are 
broad-lined SNe Ic (SNe Ic-BL). The information of the five 
LSNe is listed in Table 1. 


2.1. Modeling the Multi-band LCs of Five LSNe Using 
the Ni Model 


We first use the ?5Ni model to fit the multi-band light curves 
(LCs) of the five SNe in our sample. The details of the SN 
model can be found in Wang et al. (2023) and references 
therein. Throughout this paper, & is set to be 0.07 cm? Bm 
Definitions, units, and prior ranges of the free parameters of the 
?6Ni model are listed in Table 2. We use the Markov Chain 
Monte Carlo (MCMC) method to derive the medians and lo 
confidence regions of the parameters through the emcee 
Python package (Foreman-Mackey et al. 2013). 

The fits for the five SNe using the °Ni model are shown in 
Figure 1. The parameters and the corresponding corner plots 
are presented in Table 3 and Figures A1—A5. 

Based on the derived Mej and vpn, we can calculate the values 
of Ex of the SNe (Ex = 0.3MejVon)- We find that Ex of 
SN 2009ca, ASASSN-15mj, SN 20190md, SN 2021 bmf, and 
SN 2022ued are 1.53 x 10° erg, 1.62 x 10?! erg, 4.82 x 
10°! erg, 2.31 x 10°” erg, and 5.57 x 10?! erg, respectively. 
The derived values of Mej, Myi, Vpn, and Eg are ~ (0.7— 
8.0 Ms, ~(0.7-5.0) Mo, ~ (1.0-2.5) x 10 cm s ^!, and 
~(1.0-25.0) x 10?! erg, respectively; the mean values of these 
parameters are ~5.45 Mo, —2.79 Mo, ~ 1.73 x 10° em s^, 
and ~10.0 x 10°! erg, respectively. The ratios of Mw; to Mej 
(f) of the five SNe are —0.22-0.99. 

By comparing our derived parameters to those in the 
literature, we find that the values of derived My; and f of the 
?6Ni model are consistent with those derived from other LSNe 
using the *°Ni model, which are ~(1.0-2.0) Mo and 0.31-0.62 
(Wang et al. 2015b) or z3.2 Mẹ and ~0.31 (Gomez 
et al. 2021). 

As pointed out by Khatami & Kasen (2019) and Arnett 
(1982), Afsariardchi et al. (2021)'s model might overestimate 
the °°Ni masses of SNe. More reasonable *°Ni masses can be 
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Table 2 
Definitions, Units, and Prior Ranges of the Parameters of the 56Ni model? 
Parameter Definition Unit Prior 
56Ni modles 
Mej The ejecta mass o [0.1, 50] 
Vph The early-time photospheric velocity 10 cm s^! [0.1, 5] 
My; The °°Ni mass " [0.0, 5] 
log &4.Ni The ^-ray opacity of *°Ni-cascade-decay photons cm? g^! [-1.568, 4] 
Te The temperature floor of the photosphere K [1000, 10^] 
shift The explosion time relative to the first data day [—20, 0] 
Ahost.V The extinction in the host galaxy mag (0, 1] 
AcE The cutoff wavelength A [0, 4000] 
B The dimensionless free parameter [0, 10] 
7 /dof 
Note. 
* For SNe without UV photometric data, Acr and 8’ are set to be 3000 A and 1, respectively. 
Table 3 
Medians, 1o Bounds, and Best-fitting Values (in Parentheses) of the Parameters of the 56Ni Model 
Parameter SN 2009ca ASASSN-15mj SN 20190md SN 2021bmf SN 2022ued 
56Ni modle 
M.(Mo) 5.78015 0.84005 78*133 7.037032 6.341020 
(5.82) (0.79) (8.02) (6.35) (6.26) 
vg (10? cm s^!) 2.107003 1.991020 1.01045 2.521017 1252003 
(2.10) (1.85) (1.00) (2.47) (1.22) 
Myi(Mo) 1.267991 0.79 *015 2.267078 4.361088 4.98 +0.06 
(1.26) (0.78) (2.14) (4.77) (4.99) 
log &,, (cm g !) — 1.304983 — 1.3449 1.297133 — 1.4849 88 — 1.424533 
(—1.30) (—1.40) (—0.34) (—1.46) (—1.45) 
T(K) 6239.02::2240 2957.43* 85348 2380.15 4335-44 6173.63733438 2079.51*32041 
(6233.46) (3437.09) (1857.71) (6429.65) (2207.24) 
tsnin(day) —19.93 002 —9.95*131 —15.117}38 —9.32+963 —18.22*136 
(—19.96) (—10.18) (—15.00) (—9.36) (—18.95) 
Anos, v (mag) 0.00 5.59 0.03009 0.107045 0.374096 0.00050 
(0.00) (0.01) (0.04) (0.45) (0.00) 
Ack(À) 362528*510 3370.751324 e e e 
(3622.56) (3398.80) 
B 52208 137098 
(5.13) (1.52) ae - es 
x? /dof 22.71 1.83 0.47 2.97 8.80 
(22.69) (1.69) (0.45) (2.95) (8.65) 


Note. The values of y*/dof (reduced x”, dof = degree of freedom) are also presented. 


derived by using the equation (Khatami & Kasen 2019) 
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the mean values of 3 for Type Ib and Type Ic-BL SNe are 0.66 
and 0.56, respectively (Afsariardchi et al. 2021); the rise time tp 
and the peak luminosity L, are given by the bolometric LC 
reproduced by the best-fitting parameters of multi-band LCs 
fits. Using the equation, the ?9Ni masses of the five SNe are 
0.56 Mo, 0.44 Ms, 0.99 Mo, 2.26 Mo, and 228 Mọ, 
respectively. The values of f are ~0.10 (SN 2009ca), ~0.56 
(ASASSN-15mj), ~0.12 (SN 20190md), ~0.36 (SN 2021bmf), 
and ~0.36 (SN 2022ued), respectively. 
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Figure 1. The best fits (the solid curves) of the multi-band LCs of SN 2009ca (the top-left panel), ASASSN-15mj (the top-right panel), SN 20190md (the middle-left 
panel), SN 2021bmf (the middle-right panel), and SN 2022ued (the bottom panel) using the °°Ni model. The shaded regions indicate 1c bounds of the parameters. 


Circles represent observed data, triangles represent upper limits. 
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Table 4 

Definitions, Units, and Prior Ranges of the Parameters of the Magnetar Plus Ni, and the Fallback Plus ?9Ni Models? 
Parameter Definition Unit Prior 
Magnetar 4- 56Ni models 
Moj The ejecta mass Ms [0.1, 50] 
Po The initial period of the magnetar ms [0.8, 50] 
B The magnetic field strength of the magnetar 10^ G [0.1, 100] 
Vph The early-time photospheric velocity 10 cm s^! [0.1, 5] 
Mvi The *°Ni mass Mo [0, 0.2Mj] 
log K-, mag The ^-ray opacity of magnetar photons cm? g^! [-2. 4] 
log &4 wi The y-ray opacity of °°Ni-cascade-decay photons em? g^! [—1.568, 4] 
Tt The temperature floor of the photosphere K [1000, 10] 
shift The explosion time relative to the first data day [—20, 0] 
Ahost.V The extinction in the host galaxy mag [0, 1] 
ACE The cutoff wavelength A [0, 4000] 
B' The dimensionless free parameter [0, 10] 
x? /dof 
Fallback + °°Ni models 
Mij The ejecta mass Ms [0.1, 50] 
logL, The constant injection luminosity ergs! [48, 58] 
log tir,t The transition time day [-3, 3] 
Vph The early-time photospheric velocity 10 cm s! [0.1, 5] 
My; The *°Ni mass Mo [0, 0.2M,] 
log Ky, The y-ray opacity of fallback photons cm? g^! [-2. 4] 
log &4 wi The y-ray opacity of °°Ni-cascade-decay photons em? g^! [—1.568, 4] 
T; The temperature floor of the photosphere K [1000, 10*] 
shift The explosion time relative to the first data day [—20, 0] 
Ahost.V The extinction in the host galaxy mag [0, 1] 
ACE The cutoff wavelength A [0, 4000] 
p' The dimensionless free parameter [0, 10] 
x? /dof 
Note. 


* For SNe without UV photometric data, Acr and 3’ are set to be 3000 A and 1, respectively. 


Umeda & Nomoto (2008) suggest that the upper limit of the 
value of f of CCSNe is 0.2. Therefore, ASASSN-15mj, 
SN 2021bmf, and SN 2022ued cannot be explained by the 
°°Ni model. Although f of SN 2009ca and SN 20190md is 
smaller than 0.2, the derived ??Ni mass of SN 20190md is 
0.99 Mo, which is about twice those of some very energetic 
SNe Ic (e.g., SN 1998bw, SN 2003dh). Hence, we suggest that 
the °°Ni mass of SN 20190md is also unreasonable. Finally, the 
possibility that SN 2009ca was powered by *°Ni decay cannot 
be excluded, since its ? $Ni mass derived (0.56 Mọ) is 
comparable to those of SN 1998bw, SN 2003dh, and other 
energetic SNe Ic. 


2.2. Modeling the Multi-band LCs of Four LSNe Using 
the Magnetar Plus °°Ni Model and the Fallback Plus ??Ni 
Model 


Here, we use the magnetar plus °°Ni model and the fallback 
plus ° $Ni model to fit the multi-band LCs of ASASSN-15mj, 
SN 20190md, SN 2021bmf, and SN 2022ued. The details of the 


two models can be found in Wang & Gan (2022) and Wang 
et al. (2023) and references therein; definitions, units, and prior 
ranges of the free parameters of the two models are listed in 
Table 4. 

The fits of the two models for the four SNe are shown in 
Figure 2. The parameters and the corresponding corner plots 
are presented in Table 5 and Figures A6-A13. We find that the 
LCs of four SNe can be well fitted by the two models. The x^ / 
dof values of the best fits of ASASSN-15mj, SN 2019omd, 
SN 2021bmf, and SN 2022ued using the magnetar plus ?Ni 
model (the fallback plus *°Ni model) are respectively 0.82 
(1.10), 0.55 (0.58), 0.93 (1.00), and 1.90 (2.01), indicating that 
the magnetar plus S°Ni model is better than the fallback plus 
°°Ni model for the four LSNe. 

Based on the derived Mej and vpn, we can calculate the values 
of Ex of the SNe. For the magnetar plus *°Ni model (the 
fallback plus SON model), Ex of ASASSN-15mj, SN 20190md, 
SN 2021bmf, and SN 2022ued are 2.84 x 10°° erg (225 x 
10°° erg), 2.38 x 10°! erg (2.53 x 10°° erg), 3.35 x 10?? erg 
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Figure 2. The best fits (the solid curves) of the multi-band LCs of ASASSN-15mj, SN 2019omd, SN 2021bmf, and SN 2022ued using the magnetar plus *°Ni model (the left 
panels) and the fallback plus °°Ni model (the right panels). The dotted lines are the LCs powered by the ?5Ni, and the dashed lines are the LCs powered by the magnetar (the left 
panels) or the fallback (the right panels), respectively. The shaded regions indicate 1o bounds of the parameters. Circles represent observed data, triangles represent upper limits. 
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Medians, 1o Bounds, and Best-fitting Values (in Parentheses) of the Parameters of the Magnetar Plus 56Ni Model and the Fallback Plus *°Ni Model 


Parameter 


ASASSN-15mj 


SN 20190md 


SN 2021bmf 


SN 2022ued 


Magnetar + *°Ni Model 


M.\(Mo) 

Po(ms) 

B(10"* G) 

Vpn(10? cm s 1 
Mwi(M5) 

log Ky, ag (C0 g 5 
log yni (cm? g`’) 
TK) 

teri (day) 
Ahost,v(mag) 
er(A) 

B 

x? /dof 


Fallback + Ni model 


Me(Mc) 

logL (erg sl) 
log fy. (day) 

vg (10? em s ^J) 
Myi(Mo) 

log K». (cm? g 5 


log Ki (cm? g 


0.20033 (0.16) 
10.65 *342 (8.96) 
15.494743 (15.84) 
1.784638 (1.70) 
0.02" 9.7 (0.01) 
0.467735 (—0.53) 
1.207193 (3.11) 
3494.34 * 161539 (3890.14) 
—6.8611 78 (—6.88) 
0.29013 (0.34) 
3543.41 1037. (3598.72) 
2.507943 (2.62) 
1.59 (0.82) 


0.13065 (0.11) 
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(4.76 x 10° erg), and 1.07 x 10°° erg (6.78 x 10” erg), 
respectively. Using Equation (1) the °°Ni masses of 
ASASSN-15mj, SN 20190md, SN 2021bmf, and SN 2022ued 
derived by the magnetar plus ?9Ni model (the fallback plus °°Ni 
model) are 0.01 Mo (0.01 Mo), 0.11 M; (0.07 Mo), 0.21 Mo 
(0.01 Mo), and 0.02 Mz (0.03 Mo), respectively. 

For the magnetar plus °°Ni model, the derived values of Ma. 
Mvi, Po, B, Vpn, and Ex are ~(0.1-10) Mo, ~(0.01-0.5) Mo, 
~(2-9) ms, ~(1-15) x 10'* G, ~(0.7-2.4) x 10? cm s~', and 
~(0.10-35.0) x 10?! erg, respectively; the mean values of 
these parameters are ~3.81 Mo, ~0.18 Mo, ~5.66 ms, 
~6.18 x 10'* G, ~1.44 x 10? em s^!, and 9.07 x 10°! erg, 
respectively. 

For the fallback plus 56Ni model, the derived values of Mej, 
Mvi, Li, tito» Vpn» and Ey are ~(0.1-12) Mo, ~(0.01-0.1) Mo, 
~(0.2-8) x 10°4 erg s^!, ~(7-35) days, ~(0.7-2.6) x 10? cm 
s !, and ~(0.10-48.0) x 10°! erg, respectively; the mean 
values of these parameters are ~3.32 Mo, ~0.05 Mo, 


~4.49 x 10° erg s^, ~17.65 days, ~1.46 x IO cm s^, 
and ~1.20 x 10°* erg, respectively. 

The derived My; of the magnetar plus °°Ni model and the 
fallback plus °°Ni model are respectively ^0.01—0.21 Mo and 
~0.01-0.07 Mo, which are consistent with those of normal 
SNe Ic. Moreover, other parameters of the two models are also 
in reasonable ranges and comparable to those derived in the 
literature (e.g., Wang et al. 2015b; Gomez et al. 2022). 
Therefore, we suggest both the magnetar plus °°Ni model and 
the fallback plus ?*Ni model are reasonable. 


3. Discussion 


3.1. The Masses of the Ejecta of ASASSN-15mj and 
SN 2022ued 


The ejecta masses for ASASSN-15mj and SN 2022ued are 
0.16 Mo and 0.28 Mo (the magnetar plus °°Ni model), and 
0.1 M. and 0.23 M. (the fallback plus *°Ni model), 
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Figure 3. The Po-B distribution of the magnetar of the four LSNe (ASASSN- 
15mj, SN 20190md, SN 2021bmf, and SN 2022ued) in our sample (red dots). 
For comparison, the Po-B distributions of the magnetar of SLSNe in Nicholl 
et al. (2017) and LSNe in Gomez et al. (2022) are also plotted (blue and green 
dots, respectively). 


respectively. These values are significantly lower than the 
range of ejecta masses of SLSNe and normal Ib/c SNe in 
the literature which are 1-40 Mo (Blanchard et al. 2020) and 
1-10 Mo (Prentice et al. 2019), respectively. 

The small derived ejecta masses of the two SNe might be 
due to the absence of pre-peak data which results in poor 
constraints for the rise time and therefore the ejecta masses. 
Otherwise, the two SNe are ultra-stripped LSN like iPTF16asu 
(Whitesides et al. 2017; Wang & Gan 2022). 


3.2. Comparison with Previous Studies 


Assuming that ASASSN-15mj, SN 20190md, SN 2021bmf, 
and SN 2022ued are powered mainly by magnetars, we plot the 
distribution of Po and B (see Table 5) of the magnetars assumed 
to power the LCs of the four LSNe that we study and those of 
the magnetars associated with the LSNe studied by Gomez 
et al. (2022) and the SLSNe studied by Nicholl et al. (2017), 
see Figure 3. It can be found that the parameter space of the 
putative magnetars powering the LCs of the four SNe that we 
study overlaps both that of the LSNe studied by Gomez et al. 
(2022) and the SLSNe studied by Nicholl et al. (2017). 
Moreover, the derived upper limit of Po (~11 ms) of our 
sample is smaller than that (~23 ms) of Gomez et al. (2022), 
while the derived lower limit of Pp (~3 ms) of our sample is 
larger than that (~1 ms) of Nicholl et al. (2017). This is 
because the peak absolute magnitudes of our sample are in 
—19.5 to —20.5 mag, while the peak absolute magnitudes of 
Gomez et al. (2022) and Nicholl et al. (2017) extend to —19 
and < —21 mag, respectively. 
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Figure 4. The logL;-log ty distributions of the fallback plus *°Ni model for the 
four LSNe (ASASSN-15mj, SN 20190md, SN 2021bmf, and SN 2022ued) in 
our sample (red dots). For comparison, the logL,;-logf, distributions of the 
fallback model for SLSNe in Moriya et al. (20182) are also plotted (blue dots). 


Assuming that ASASSN-15mj, SN 20190md, SN 2021bmf, 
and SN 2022ued are powered mainly by fallback, we plot the 
distribution of L4, and ty (see Table 5) of the four LSNe and 
those of SLSNe studied by Moriya et al. (20182), see Figure 4. 
It can be found that the parameters of the four SNe we study are 
located in the lower right of the parameter space (larger ty and 
lower Lı) while the parameters of SLSNe studied by Moriya 
et al. (2018a) are located mainly in the upper left of the 
parameter space (smaller f and higher L,). It should be noted 
that, however, there is a small overlap between our sample and 
the sample of Moriya et al. (20182). 


3.3. The Accretion Masses and the Validity of the 
Fallback Plus *°Ni Model 


For the fallback plus S°Ni model, the accretion mass Macc 
can be written as (Moriya et al. 2018a) 


Mace? = 2:51 5 (2) 


here, 7) is the efficiency of converting accretion to input energy 
which can be from ~0.001 (Dexter & Kasen 2013) to ~0.1 
(e.g., McKinney 2005; Kumar et al. 2008; Gilkis et al. 2016). 
Using the equation, the accretion masses of the four LSNe are 
0.0005—0.0547 Mo, 0.0043-0.4345 Mo, 0.0043-0.4378 Mo, 
0.0084—0.8472 Mọ, respectively. The accretion masses are 
significantly smaller than the derived *°Ni masses if 7 is set to 
be the upper limit (0.1). In contrast, the accretion masses are 
larger than the derived ?9Ni masses if 1] is set to be the lower 
limit (0.001). Therefore, the validity of the fallback plus 56Ni 
model depends on the value of ņ and favors a large 7) value. 
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4. Conclusions 


In this paper, we study five LSNe with Mpeak zz —19.5 to 
—21mag by fitting their multi-band LCs with different energy 
source models. We find that the °°Ni model cannot account for 
the LCs of four SNe (ASASSN-15mj, SN 20190md, 
SN 2002ued, and SN 2021bmf) in our sample, since f of the 
three SNe (ASASSN-15mj, SN 2002ued, and SN 2021bmf) are 
larger than the upper limit (0.2) of f value of CCSNe, while the 
derived °°Ni mass of SN 20190md is 0.99 Mo, which is about 
twice those of some very energetic SNe Ic. For comparison, the 
derived ?9Ni mass (0.56 Mo) of SN 2009ca is comparable to 
those of energetic SNe Ic, indicating that it might be powered 
by the ?9Ni cascade decay. This indicates that some LSNe 
might be powered by ?9Ni cascade decay, though most of them 
need additional energy sources. 

We use the magnetar plus °°Ni model and the fallback plus 
°°Ni model to fit the LCs of the four LSNe that cannot be 
explained by the ?9Ni model, finding that the two models can 
account for the four SNe. The X? /dof values of the best fits of 
the two models suggest that the magnetar plus S°Ni model is 
better than the fallback plus °°Ni model for the four LSNe. 

For the magnetar plus S°Ni model, Mej, Mwi, Po, B, Vpn, and 
Ex are ~(0.2-10) Mo, <0.3 Mo, ~(2-9) ms, ~(1-15) x 10^ G, 
~(0.7-2.4) x 10? cm s~!, and ~(0.10-35.0) x 10°! erg, respec- 
tively. These values are comparable to the typical values of the 
magnetars supposed to power the LCs of LSNe (see, e.g., Wang 
et al. 2015b; Gomez et al. 2022) and some SLSNe (see, e.g., 
Nicholl et al. 2017). For the fallback plus °°Ni model, Mej, Mw, 
Li, fijfo; Vpn, and Ex are ~(0.1-12) Mo, ~(0.01-0.1) Mo, 
~(0.2-8) x 109^ erg s^!, ~(7-35) day, ~(0.7-2.6) x 10? cm 
s^ !, and ~(0.10-48.0) x 10°! erg, respectively. 

It should be noted that, however, the derived Mej of 
ASASSN-15mj and SN 2022ued are respectively 0.16 Mo 
(0.11 Mo) and 0.28 Mo (0.23 Mo) for the magnetar plus *°Ni 
model (the fallback plus ?9Ni model), which are comparable to 
those of ultra-stripped SNe. This might be due to the absence of 
pre-peak data which prevents us from better constraining the 
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rise time and therefore the ejecta masses; or the two SNe are 
ultra-stripped LSN like iPTF16asu. 

Finally, we find that the distribution of Po and B of the 
putative magnetars powering the LCs of the four SNe overlaps 
those of the magnetars of LSNe and SLSNe in the literature, if 
they are (mainly) powered by magnetars. Similarly, the 
distribution of Lı and t of the fallback powering the LCs of 
the four SNe overlaps those of the fallback of SLSNe in the 
literature, if they are (mainly) powered by fallback. 

We suggest that the magnetar plus ??Ni model is more 
reasonable than the fallback plus °°Ni model, since the validity 
of the fallback plus °°Ni model depends on the value of 7 and 
favors a large 7) value, and the magnetar plus °°Ni model yields 
smaller x /dof values. It should be pointed out that, however, 
the fallback plus *°Ni model is still a promising model that can 
account for the four SNe in our sample as well as other LSNe. 
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Appendix 


Figures Al—A5 are the corner plots of the °°Ni model of 
N 2009ca, ASASSN-15mj, SN 2019omd, SN 2021bmf, and 
N 2022ued, respectively. Figures A6—A9 are the corner plots 
f the magnetar plus °CNi model of ASASSN-15mj, 
N 20190md, SN 2021bmf, and SN 2022ued, respectively. 
igures A10-A13 are the corner plots of the fallback plus 
?5Ni model of ASASSN-15mj, SN 2019omd, SN 2021bmf, and 
SN 2022ued, respectively. 
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Figure A1. The corner plot of the *°Ni model for the multi-band LCs of SN 2009ca. The solid vertical lines represent the best-fitting parameters, while the dashed 


vertical lines represent the medians and the 1o bounds of the parameters. 
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Figure A2. The corner plot of the *°Ni model for the multi-band LCs of ASASSN-15mj. The solid vertical lines represent the best-fitting parameters, while the dashed 


vertical lines represent the medians and the 1e bounds of the parameters. 
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Figure A3. The corner plot of the *°Ni model for the multi-band LCs of SN 20190md. The solid vertical lines represent the best-fitting parameters, while the dashed 


vertical lines represent the medians and the 1o bounds of the parameters. 
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Figure A7. The corner plot of the magnetar plus *°Ni model for the multi-band LCs of SN 20190md. The solid vertical lines represent the best-fitting parameters, 
while the dashed vertical lines represent the medians and the 1c bounds of the parameters. 
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Figure AS. The corner plot of the magnetar plus *°Ni model for the multi-band LCs of SN 2021bmf. The solid vertical lines represent the best-fitting parameters, while 
the dashed vertical lines represent the medians and the 1c bounds of the parameters. 
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Figure A9. The corner plot of the magnetar plus °°Ni model for the multi-band LCs of SN 2022ued. The solid vertical lines represent the best-fitting parameters, while 


the dashed vertical lines represent the medians and the 1o bounds of the parameters. 
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Figure A10. The corner plot of the fallback plus °°Ni model for the multi-band LCs of ASASSN-15mj. The solid vertical lines represent the best-fitting parameters, 
while the dashed vertical lines represent the medians and the 1c bounds of the parameters. 
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Figure A11. The corner plot of the fallback plus *°Ni model for the multi-band LCs of SN 20190md. The solid vertical lines represent the best-fitting parameters, 
while the dashed vertical lines represent the medians and the 1c bounds of the parameters. 


20 


Research in Astronomy and Astrophysics, 23:125002 (23pp), 2023 December 


Maj (Mo) = 14. 
T 


Igu (ergs?) 


+6.84 
367550 


gu (ergs ?) = 54.4709 


og(T«)(days) 


log(Tir)(days) 


e) = 0.084932 


Bai et al. 


Anost (mag) = 0.697338 


ri 
= 
z 
zo 
log(Ky faliback) = 1.294385 
3 
ud 
> 
É 
o 
2 
° 
log(ky ni) = 0.704224 
CE 
= 
EÀ 
D 
2 
s 
T: (K) = 6868.95*59632 
S 
e 
gS | 
E 
e 
= 0.66. 
= 9.907088 
i 
E 
= 
o 
2, 
5 
a 
v 
S 
~ o? 
D 
8 
E 
iy 
bd 


Maj (Mo) 


lgu (ergs) 


log(Tr)(days) 


log(ky, fallback) 


log(ky, ni) 


Ti (K) 


tonire(days) 


nost (mag) 


Figure A12. The corner plot of the fallback plus 56Ni model for the multi-band LCs of SN 2021bmf. The solid vertical lines represent the best-fitting parameters, while 
the dashed vertical lines represent the medians and the 1o bounds of the parameters. 
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Figure A13. The corner plot of the fallback plus *°Ni model for the multi-band LCs of SN 2022ued. The solid vertical lines represent the best-fitting parameters, while 
the dashed vertical lines represent the medians and the 1c bounds of the parameters. 
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